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The n-dimensional hypercube quantum random walk (QRW) is a particularity appealing example 
of a quantum walk because it has a natural implementation on a register on n qubits. However, 
any real implementation will encounter decoherence effects due to interactions with uncontrollable 
degrees of freedom. We present a complete characterization of the mixing properties of the hypercube 
QRW under a physically relevant Markovian decoherence model. In the local decoherence model 
considered the non-unitary dynamics are modeled as a sum of projections on individual qubits to 
an arbitrary direction on the Bloch sphere. We prove that there is always classical (asymptotic) 
mixing in this model and specify the conditions under which instantaneous mixing always exists. 
And we show that the latter mixing property, as well as the classical mixing time, depend heavily 
on the exact environmental interaction and its strength. Therefore, algorithmic applications of the 
QRW on the hypercube, if they intend to employ mixing properties, need to consider both the walk 
dynamics and the precise decoherence model. 



I. INTRODUCTION, PREVIOUS WORK, AND OUR WORK 

Quantum walks [I] play a prominent role in the design of quantum algorithms. Their distinction 
from classical random walks lies in their potentially faster mixing and hitting times. The under- 
lying dynamics of quantum walks can be either continuous or discrete, and even though these two 
representations have some properties in common, analysis usually demands different techniques and 
tools. Despite their different origins, discrete and continuous time quantum walks can be precisely 
related to each another |42j . 

The theoretical properties of quantum walks on general graphs have been outlined in Ref. |10j . 
and the remarkable result by Szegedy [9J explains a general framework for the proper quantization 
of any Markov chain algorithm. The dynamics of quantum walks have been analyzed, for example, 
on the line [TT], circle [IS], hyperlattice [T3] and hypercube [2]. Shenvi et. al. [7] proved that the 
discrete quantum walk can be used in a search algorithm and Ambainis used it for the best known 
algorithm for element distinctness [5]. The continuous-time quantum walk was originally proposed 
by Childs, Farhi and Guttman [11 [5] as an algorithmic primitive. Childs et. al. used a quantum walk 
to prove the best known results for the separation between quantum and classical query complexity 
[B], and a highly efficient algorithm for NAND formula evaluation [35] . 

Central to the algorithmic application of both classical and quantum random walks are their 
mixing characteristics. For a classical random walk that has a unique steady state, the mixing 
time characterizes the convergence of the walk to this steady state. In the quantum case, unitarity 
prevents the walk from reaching a steady state. This has led to alternative notions of mixing for 
quantum walks. One is to define an instantaneous mixing time, as the first instant the probability 
distribution of the walker's location on the graph is e-close to the uniform distribution. Another 
sensible definition for the mixing time of a continuous quantum walk, although one that is dependent 
on the initial state, is based on a limiting value of a time-averaged probability distribution |10j . 

Quantum systems are very susceptible to imperfections and interactions with their environment, 
both of which cause decoherence. Sufficient decoherence can remove any potential benefits from the 
quantum dynamics. Investigations to date have either used Markovian models for the environment, 
with the environment monitoring the walker position or state of the 'coin' driving the walk |16l 1171 
HHJ [191 EOl EQl EI] ; or imperfect evolution [2j E3 Ell E3 EE] [30] , such as broken edges. 
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While decoherence is nominally the nemesis of quantum information processing, it has been argued 
that decoherence can in fact be 'useful' in the context of quantum walks |16j . Decoherence can be 
used to force a quantum walk to mix to a uniform distribution, and in this paper we will illustrate 
this for the quantum walk on the hypercube. Similar results have been shown for quantum walks 
on the line [H] and TV-cycle [21] using weak measurements of the walker's position. Additionally, 
Richter has provided elegant analyses of the general phenomena of quantum speedup of classical 
mixing processes under a restricted decoherence model |29j . 

In ref. |30j . Marquezino et al. examined the discrete-time quantum walk on the hypercube, and 
derived the limiting time-averaged distribution in the coherent case (no decoherence). The mixing 
behavior, both to this distribution and the uniform distribution, were considered for a coherent walk 
as well as under the decohering effects of randomly breaking links in the network. In the decoherent 
case, the walk was shown to approach the uniform distribution. Interestingly there is an optimal 
decoherence rate which provides the fastest convergence. A similar effect was found for the iV-cycle 
in |16j . Below we show that this can also be true for the continuous-time quantum walk on the 
hypercube, but not always. 

Hitting times 1 and instantaneous mixing times for the continuous-time version of the hypercube 
quantum walk with decoherence were calculated recently by Alagic and Russell |21j . Analytical 
results were derived by exploiting the representation of the quantum walk on the hypercube as a 
set of non-interacting qubits; a simple example of how spin networks may be mapped to quantum 
walks [31] ■ We extend these results and provide a complete characterization of the mixing time in 
terms of the decoherence intensity and projection direction. In our case the projection direction 
of the decoherence operator can be in any arbitrary direction in contrast to the previous analysis 
which was restricted to decoherence in the computational basis {|0), |1)}. Moreover, we show that 
randomizing the direction of decoherence leads to the depolarization channel that, regardless of the 
decoherence intensity, shows universal mixing behavior. 

The paper is organized as follows: we first describe the quantum walk on the hypercube and 
its mapping to a register of qubits in section [IT] Then in section |III| we introduce our decoherence 
model and explicitly show that it is the continuous-time analogue of the standard discrete-time 
projection model. The definitions of various mixing times are introduced formally in section |IV| 
With our quantum walk and decoherence model the hypercube quantum walk is separable as a 
product over single qubits, and we need only consider single qubit dynamics, which we discuss in 
section [V] Analysis of simple channels, randomized projections and numerical results are presented 
as subsections of section |V] 



Continuous-time quantum walks [1] are defined over an undirected graph with N = 2™ nodes, each 
labelled by an integer i € [0,7V — 1]. These walks can in general be described by the Hamiltonian 



where each node i corresponds to the quantum state, \i) = cj 1 0) , and [ij] denotes connected nodes 
i and j. The state \i) thus corresponds to a 'particle' located at node i. The first term in is 
a 'hopping' term with amplitude Ajj(t) between nodes i and j; the second describes 'on-site' node 
energies tj (t) . Both these terms can depend on time. One can simplify (JIJ by dropping all the on-site 
energies, and by making all the intcrnodc hopping matrix elements the same, i.e., A$j — > A,V{i, j}. 

The structure of a hypercube is particularly appealing; it is an n— regular graph that is the 
underlying model for many computational problems. The nodes of this graph can be represented as 
basis vectors \v) £ {|1), |0)}® ra in C 2 dimensional Hilbert space - the same Hilbert space describing 



The hitting time is defined as the first time a given vertex, or set of such, is reached. 



II. QUANTUM WALK ON THE HYPERCUBE 





(1) 
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n qubits. Each node is labelled by a binary string representing a multi-qubit state, i.e. \z) = 
\z\Zi . . . z„) = \zi) ® l^) ® ... <8> |z n ), where the z^s are or 1; pairs of nodes with a Hamming 
distance of 1 (the number of bits that must be flipped to obtain one from the other) are connected to 
give the hypercube. In this way n qubits describe a quantum walk over a TV-dimensional hypercube, 
which takes place in information space. This is described by the simple qubit Hamiltonian 

n 

ff = A^>i, (2) 

<=i 

representing a set of non-interacting qubits, each evolving under the Pauli a x = ( ^ q J operator. 

Note that we do not scale the rate of Hamiltonian evolution with the number of qubits (e.g. H oc A/n 
for some A independent of n) as was done in some prior work (e.g. Ref. [14]). We believe this scaling 
is artificial and unnecessary for a quantum walk implemented using a qubit register. 

In this case the unitary dynamics is trivially solvable; for the walker initialized at the z = corner, 
the probability of being at some site z is 

n 

P,(t) ^ (z\g(t)\z) = (z\ (gje-^lOMOle^ 4 \z) 

i=l 
n 

= {A 0[cos(At)|O)+sin(At)|l))(cos(At)(O|+sin(At)(l|] \z) 

i=l 

= cos 2 ™ (At) sin 2 ™ 1 (At) 

(3) 

where no is the number of O's, and ri\ the number of l's appearing in z, and g(t) = e~ tHt \z — 0) (z = 
0|e lff * is the density matrix of the qubit register at time t. The state of the qubit register is always 
pure, but we specify it with a density matrix to be consistent with what follows. 

This mapping to a qubit model is not only useful conceptually, but it is highly suggestive of a 
potential physical implementation of the hypercube quantum walk. If the system is not closed and 
it is exposed to measurement or an environment that it interacts with, the quantum dynamics are 
in general far more complicated and we will examine this now. 



III. DECOHERENCE MODEL 



To study decoherence of quantum walks it is common to begin with a discrete-time quantum walk 
and model the decoherence as a sequence of weak measurements on the walker [T5J [TBI HZ1 UM ■ When 
the result of the measurement is ignored (i.e. lost to the environment), the nonunitary process is 
described by: 

g n+1 = (\-p)U en tf + P Y J M a {U Qn tf)Ml (4) 

a 

where the measurement is given by the POVM {M a }, such that ^ Q M a = 1 and M a > 0, occurring 
with probability p at each time step; U describes the unitary evolution of the quantum walk. This 
model is equivalent to a memoryless environment, unperturbed by the system. We will consider the 
continuous-time analogue of this process, which we derive below. 

Claim: The discrete-time weak measurement model of the system dynamics, 

Qt+r = (1 - yr)U r Qtl4 + rr^2M a [U r gtl4]Mi . ( 5 ) 

a 

is equivalent, in the limit r — ► to the master equation 



g(t) = -i[H,g{t)]+'r^2v[M a ]Q(t) 

a 



(G) 
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if the measurement rate 7 is such that for a time-step of duration r, p = 7T, unitary evolution 
U T = exp(-iHr). Here the superoperator T>[X]g = Xg{t)X^ - \ (X^Xg{t) + g(t)X*X) for any 
operator X. 

Proof: For small r we expand the exponential to first-order such that Eq. §5§ becomes 

g t+T « (l-~,T)(l-iTH) gt (l + iTH)+jTj2M a (l-iTH)Q t (l+iTH)Ml (7) 

a 

Now, dropping second order terms in the parameter r, 

ft+r ~ ft + irfti? - irHg t ~ jrg t + 7T ^M aQt Ml (8) 

CK 

= g t + iT Qt H - ir/fft +jrJ2 Ul a g t Ml - ^M\,M a Q t - J . (9) 

Dividing by r we obtain, 

gt+T r ~ gt = ft] + 7 E ( M ^* M 1 " 5 WlM a g t - g t MlM a ] \ , (10) 

then taking the limit r — > 0, 

g(t) = -i[H, g(t)} + 7 ^ f M a g(t)A4 - |[ikfj M a(? (t) - e(t)MjM a ] ) (11) 

= -i[fr, e (t)] + 7 5^i>[M a ]e(t) , (12) 

a 

as claimed. ■ 

The interaction between the system and the environment can be represented in various ways, 
through the choice of the set of {M a }. In the quantum walk literature, the standard choice is 
the walker location, i.e. M a — = P,*, the projectors onto the graph node states. However, 
when the quantum walk on the hypercube is implemented via a set of qubits, as we describe here, 
position measurements corresponds to a computational basis measurement of the state of every qubit 
simultaneously. This implies a physically unrealistic, multi-qubit measurement/interaction with the 
environment. 

If the quantum walk were to be implemented using a qubit register, a more physically realistic 
decoherence process is described by single-qubit projective measurements. The obvious choice, as 
an analogue to the location measurements, are measurements that are projections onto single qubit 
computational basis states; this is the choice considered in [3T]. We can generalize this to single- 
qubit projective measurements onto arbitrary antipodal points on the Bloch sphere. We express this 
as: 

I + ~r ■ a _ I + (sin 9 cos tp) a x + (sin 9 sin tp)a y + (cos 9)a z „2x2 
F oC r J = o = 2 

p x (^) = 1 ~ T. ' ° - 1 ~ (siQ^cosy)^ - (sin g sin y)^ - (cos 9) a z c C 2x2 ^ 

where I is the two-dimensional identity matrix, ~r — (r x ,r y ,r z ) — (sin 9 cosy, sin 9 sin tp, cos 9), for 
9 G [0,7r], tp £ [0,27r], defines what we refer to as the "decoherence axis" or measurement projection 
direction. By a = (a x , o~ y , a z ) we denote the three dimensional matrix vector composed of the three 
nontrivial Pauli matrices. This decoherence model is a generalization of what has been referred to 
as the subspace projection decoherence model [3]. 

We consider the continuous quantum walk described by the master equation Eq. ([6]), with a 
Hamiltonian given by Eq. ^ and measurement operators given by 

M* = I®...®P a ("r)®...®I 5 (14) 

where a = {0, 1} and the projector is on the k th qubit. Note that ^2 a = I® n . The qubit register 
evolution equation can then be written as: 

n n 1 

m = e -iA [at g (t)} + 7 E E 13 M <?(*)• ( 15 ) 

k=l fc=lQ=0 
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This can alternatively be written as: 



e(t) = K g(t)] + IY. V ^-^] ( 16 ) 

fc=i fc=i 



where the sum is over qubits. None of the qubits are interacting and therefore this master equation 
has a separability property that allows us to treat the dynamics as n single qubit density matrices 
undergoing the evolution: 

p k {t) = -iA [a x ,p k (t)] + lv\t- a] Pk (t), (17) 

and g — Pk- See the appendix for an explicit derivation of this. This property allows one to 

analyze the full system dynamics by looking at single qubit dynamics, since the dynamics of the 
system is just the tensor product of the dynamics of individual, non-interacting qubits. Furthermore, 
the evolution equation for pk is the same for all k, and so we will drop the subscript when referring 
to single qubit dynamics. 

We will investigate how changes in the single qubit dynamics affect properties of the quantum 
walk on the hypercube. After formally defining the mixing time of a quantum walk in the next 
section, we provide a complete characterization of the mixing time in terms of the quantum channel 



in Eq. (17), i.e. dependence on the physical rates and the direction of the decoherence axis, r 
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IV. MIXING TIME 



To identify the physical quantities of interest, we return to a principal motivation for considering 
random walks (quantum or classical). In computer science, the most efficient solution to many 
problems is given by a probabilistic algorithm, where the correct answer is attained with high- 
probability if the space of solutions is sampled with a well-chosen sampling distribution. Generating 
the correct sampling distribution is often a matter of mapping the uniform distribution into the 
desired one, and therefore generating a truly uniform distribution is an important problem. 

There arc several different definitions in the literature that have been used as a measure of mixing 
time. For completeness will briefly list them all here. For our purposes the notions of instantaneous 
mixing and classical mixing will be the most important ones. 

Instantaneous mixing is defined as the first time instant at which the probability distribution of 
the walker's position is sufficiently close to uniform distribution: 

M insttE = min{t | \\P(x,t) - P u \\ tv < e} (18) 

where P(x, t) is the probability of obtaining clement x € X (X is the space of events we are 
sampling from, which in the case of random walks is the space of the walker location parameter) 
at time t, and P u is the uniform distribution over X. \\ ■ \\ tv is the total variation distance over 
probability distributions (we will restrict our attention to finite sample spaces). This definition is 
mostly used in idealized continuous quantum random walks where no decoherence effects are present. 
Although formally present by this definition, mixing in continuous quantum random walks without 
decoherence is only an instantaneous phenomenon. The ability to harness this instantaneous mixing 
is still questionable. 

Average mixing is based on the time-averaged probability distribution, that even for unitary 
quantum walks is shown to converge jTUj. In the continuous-time case the time-averaged probability 
distribution for the state x is defined as: 

P(x,t) = i / P(x,t)dt, (19) 
t Jo 

We can define the corresponding time-averaged mixing as: 

M avg , e = min{T | Vr > T : \\P(x, r) - P tt || to < e}. (20) 



2 Note that varying r is equivalent to changing the basis in which node states are encoded while the decoherence 
axis is kept fixed. 
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This time-averaged distribution can be sampled from by selecting t uniformly in [0,T], running the 
quantum walk for time t, and measuring the walker position. 

Classical (asymptotic) mixing is the quantity originally used in classical random walks and it 
defines mixing as the time after which the register's distribution is desirably close to uniform: 



M class , £ = min{T | Vt > T : \\P(x,t) - P u \\ tv < e} 



(21) 



This definition characterizes the time it takes for the probability of finding the walker at a particular 
location to be distributed uniformly across the entire sample space, X. 

The mixing time is a well defined quantity for classical random walks because there exists a 
stationary distribution for classical random walks over any connected, non-bipartite graph [37], 
however for continuous quantum walks this is not necessarily the case; unitary dynamics means the 
probability distribution over the graph nodes oscillates, and therefore never converges to the uniform 
distribution. We shall return to this issue below. But first, let us examine the total variational 
distance in the context of a quantum walk on a hypercube. 

For a hypercube quantum walk implemented using qubits, the walker location is encoded into the 
value of the qubit register in the computational basis. Therefore, the sample space in this case is 
the space of binary strings of length n, and the probability of measuring any register value (walker 
"location" ), z is: 



P(x,t) = (l-po(t)) k po(t) n 



-k 



(22) 



where k is the Hamming weight of the binary string z (i.e. number of ones in z), and Po(t) is the 
probability of a qubit value being at time t. Note that: 



Po(t) 



1 
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where (A (t)) = tr(Ap(t)) for any operator A. The total variational distance in this case is: 



\\P(z,t) - P u \\ tv 



E 



i 



2™ 

Po(t)) k Po(t) n - k 



(23) 



We can bound this variational distance using the Hellinger distance on distributions |21j . The 
Hcllinger distance, for two distributions V(x) and Q(x) (both defined over the same sample space 
X), is defined as: 



(24) 



Its usefulness for us comes from its relation to the total variational distance: 

\\V - Q\\ tv < 2H(T, Q) <2\\V- Q\\\i 2 

Therefore in our case, 

\\P(z,t) - P u \\ 2 tv < 4H(P(z,t),P u ) 2 



= 4- 4 Wp(^=4- 4 W? 



4-4 




(25) 



(l-Po(t))*po(*) 



n—k 



1(^/1 _( <72 ( i )) + v /i + ( <72 ( i ))) n 



(26) 



Note that this is a positive quantity because we only take the positive branch of the square roots, and 
it attains its minimum value of zero when (a z (t)) = 0. Furthermore, from this bound we see that the 
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variational distance for the distribution of the register of n qubits (from the uniform distribution) 
is small exactly when the variational distance for the distribution of a single qubit is small. In fact, 
when (a z ) is small, we can expand the square roots to second order in this quantity to get: 

\\P(z,t) - P u \\l < 4 - 4 (l - (27) 



Given this fact, we will concentrate on the distribution for a single qubit in the following, and appeal 
to the fact that the variational distance for the entire register from the uniform distribution is small 
precisely when (a z (t)) 2 is small for a single qubit. 

Finally, although our intention in this work is not to characterize the mixing behavior of the 
random walk as a function of n, the number of qubits, we note that the above inequality allows us 
to easily place bounds on this behavior. Consider the time at which the total variational distance 
of the register of qubits is at most e; the above bound tells us that this is when: 

„ln(l-fef>!)=l„(l-|) (28, 

Expanding both logs in the small parameters e and a z (t)} 2 /8 (because (cr z (t)} will be small whenever 
this equality is fulfilled) we get (a z (t)) = y^f • Now, we shall show in the following sections that 

the behavior of the single qubit (<J z (t)) as a function of time always has an exponentially decaying 
envelope when there is non-zero decoherence (7 7^ OV Therefore, to bound the time for asymptotic 



mixing we can substitute (<J z (t)} — Ce~ rt into Eq. (28 1, and get the following bound on the behavior 
of register: t m j X = O(lnn). This scaling behavior has been noted elsewhere and will not be the 
subject of study in this work. Instead, we will investigate how the mixing of a single qubit is 
affected by the physical parameters in the problem, particularly the decoherence parameters. In 
effect, we will study how the values of constants C and r in the exponential envelope above are 
determined by the physical parameters. 



V. SINGLE-QUBIT DYNAMICS 

In order to solve for the single qubit dynamics we use the following parametrization of the single 
qubit density operator, 

„ m _ l ( 1 + M*)> Mt)) - i{a y {t)) \ , , 

P{t) ~2\{cr x {t))+i{a y {t)) l-(a z (t)) ) ' ( ' 

This operator is completely described by the Bloch vector (<r) = ((<r x ), {<r y ), (c z )) = (x 7 y,z). We 
can then derive the matrix equation for the Bloch vector of a single qubit when the evolution is 



given by Eq. (17), repeated here for clarity: 



p(t) = -iA [a x ,p(t)} - lp{t) + \ (r • a) p(t) (r ■ ^ . (30) 

Multiplying this equation by each of the Pauli matrices and taking the trace we obtain the 3x3 
parametrized linear system equation: 

(*)) = |tr(^p(i)) - tr(^p(i)) - A<7? (t)) , (31) 

where, 

/ l( r l - !) jr x r y jr x r z \ 

A = jr x r y 7(7-2 _ y 1Tyrz _ 2A . (32) 
\ ir x r z jr y r z +2A ^{r\ - 1) / 

The solution of this system is (It (t)) = exp(Af)(7r (0)), where ( a (0)) is the Bloch vector of the 
initial state |0). The dynamics of the qubit, and in turn the entire register implementing the quantum 
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walk, is completely determined by the properties the matrix A. We now examine the key properties 
of this matrix. 

Firstly, by the Routh-Hurwitz criterion [57], the matrix A has eigenvalues that lie in the left half 
of the complex plane for positive 7, A and all decoherence axes except for when r x = 1. This singular 
case represents a channel where the Hamiltonian dynamics and the decoherence dynamics commute. 
This case is easy to solve for explicitly (we present the solution in 
fairly uninteresting. For all other parameter regimes, the Routh-Hi 



V A I and the dynamics for it are 



urwitz criterion tells us that Eq. 



(31 1 is a strictly stable system. 

The eigenvalues of A can be determined from its characteristic equation: 

A 3 + 2 7 A 2 + ( 7 2 + 4A 2 )A + 4 7 A 2 (1 - r 2 ) = (33) 

Interestingly the eigenvalues depend only on r x , 7, A. The most convenient form for the solutions to 
this cubic equation can be found by mapping the equation into a third order Chebyshev polynomial 
and using the Chebyshev cube root. Then the solution can be written in closed form as: 

„ /V - 12A 2 1 . NN 2 7 

Ai = 2W cos(- arccos(m)) - — 

Pf _ 12 A2 l 2 7 
A 2 = -2y cos(- arccos(-m)) - — 

A3 = -A1-A2-27, (34) 

where 

m = ( 7 2 _ i 7 2 A 2 )3 / 2 ^ + 18A 2 (3r 2 - 1))). (35) 

These eigenvalues fall into one of two classes, depending on the values of 7, A, and r : one real and 
two imaginary eigenvalues, or three (with possible repetitions) real eigenvalues. 

A is not a symmetric matrix and is therefore not generally diagonalizable. A could be not 
diagonalizable if it has repeated eigenvalues (i.e. a degenerate eigenspace). We will see below when 
we perform a more detailed analysis of the eigenvalues of A that this only occurs in a vanishingly 
small parameter range which will not be of interest to us. Therefore we will effectively treat A as 
diagonalizable. 

Given these properties of the matrix A, let us return to the mixing behavior of a single qubit 



under the dynamics given by Eq. (31 ). Assuming that the initial state of the whole register of the 
system is in the state |0)®" at t = the (er z (t)) component of each individual qubit is: 

(a z (t)) = (0 l)e A V (0)) = (0 l)e At (0 1) T = (e At ) 33 (36) 
For diagonalizable A, with eigenvalues Xj, e At can be written as: 



Using this expansion, 



M*)>= E 7x T^Tx r ^[-4A 2 +7 2 (l-r 2 ) + (A 7r2 +A. 3 ) 7 (l-r 2 ) + A T2 A^] 

ir(Ai,A 2 ,A 3 ) ' 1 _ A 7 r 2 )(A 7 r i — ^ 3 ) 

(38) 

where 7r(Ai Aa,.^) denotes the three cyclic permutations of (Ai,A2,A3). From this expression, we 
can see that (a z (t)) depends on all the free parameters in the system: r x ,r z ,^, A (even though the 
eigenvalues only depend on r x ,j and A). 

This expression for (a z )(t) tells us something crucial about the QRW. The exponential envelopes 
e A,r i * , and the fact that A is a strictly stable matrix (has eigenvalues in the left half of the complex 



plane), imply that (a z ) ^°°> 0. Hence, by Eq. (27), the limiting distribution for the quantum walk 



for all r and 7 > is the uniform distribution. Thus, the decoherence ensures that the random 
walk mixes to uniform given sufficient time. 
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Finally, we note that although the above convergence argument was for the particular initial state, 
|0)® ra , the strict stability of the dynamical matrix A (for all parameters except r x = 1) allows us to 

state more generally that there is an initial state independent steady state (~a) *^°°> 0, which of 

course implies (a z ) t ^°°> 0. 



A. Special Cases: Simple Channels 



In this section, to illustrate the utility of our approach we exactly solve the dynamics for several 
simple single-qubit channels. 

For r x = 0, the dynamics is described by the matrix 




A = 7 (r^ - 1) iry r z - 2A . (39) 



So (<j x (t)) = e -7 (<7 a (0)), and we have a pair of coupled linear differential equations describing 
the motion in the y-z plane. The motion in this plane is analogous to a damped simple harmonic 
oscillator with natural frequency 2A and damping rate 7; for 7 < 4A the system is underdamped, 
resulting in decaying oscillations around the origin, while in the overdamped regime, 7 > 4A, we 
see exponential decay. In analyzing the mixing properties of the hypercube quantum walk, we 
need only consider the behavior of the z-component of the Bloch vector. When the initial state is 
(cr (0)) = (0, 0, 1) T , in the underdamped case, 7 < 4A, we have, 

(a z (t)) = e~ 7t/2 ( cos u>t + 7 ^ ~ ^ smut) . (40) 



where w = y/\j 2 - 16A 2 |/2. In the overdamped regime, 7 > 4 A, 

M*)> = , 2 * A2 ((*+ + irl)e X+t - (X- + irlV-t) , (41) 
V7 2 — \bA z 

where A± = (—7 ± \/"/ 2 — 16A 2 )/2. Finally, in the critically damped case 7 = 4A, [a z {t)) = 
(2A(2r2 - l)t+ 1) e- 2At . 

This result applies to two standard, single-qubit decoherence channels [33]; the phase- flip, or 
dephasing, channel where ~r — (0,0, 1) and the bit-phase-flip channel r = (0, 1,0). 

In the opposite case, where r = (1,0,0) (the bit-flip channel), 




A = — 7 -2A , (42) 



and the solution for z(t) ,with initial condition (~a (0)) = (0, 0, 1) T , is simply 

(a z (t)) = e" 7 ' cos 2 At. (43) 
Finally we examine the depolarizing channel |33j . which although is not an instance of the class of 



channels described by ( 17 1, can be obtained by a randomization of such channels (see the Appendix 
for the detailed derivation of this). The depolarization channel corresponds to a randomization of 
the decoherence axis, r , and is described by a dynamical matrix: 







.21 

A = I 3 -2A I . (44) 



3 

2A 



The solution for the z-component is: 



(cr 2 (i)) = e"^ t cos2A£. 



(45) 
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B. The Complete Classification of Mixing Behavior 

For dynamics under a general single qubit channel, the convergence of the (a z (<)) to the limiting 
value can be essentially of two different types, depending on the eigenvalues of A: exponential 
decay or dampened oscillator decay. We now characterize the type of decay in terms of the physical 
parameters of the qubit model, and then use the decay types to make conclusions about the mixing 
properties of the quantum walk. 

As we have already seen the characteristic equation for the matrix A is: 



A 3 + 2 7 A 2 + ( 7 2 + 4A 2 )A + ^A 2 t] = 0, 



(46) 



denoting r\ 



For a cubic equation with real coefficients the solutions are either: (1) three 



real, possibly repeated, roots, or (2) one real root and two complex conjugate roots. According to 
Eq. (38) these two classes of roots give rise to two different types of convergence of (a z )(t) to zero: 
exponential decay or dampened oscillator decay. To determine the parameter dependence of these 
two types of convergence behavior, we can use the discriminant of (46 1: 



,2\2\ 



(47) 



A(rj) = 432 7 2 AV - (16 7 4 A 2 + 576 7 2 A 4 )r/ + (4( 7 2 + 4A 2 ) 3 - 4 7 2 ( 7 2 + 4A 
There are two fundamentally distinct parameter regions: 

(i) (Zeno region) A(rj) < =>• 3 real roots, that are distinct unless, 
A(rj) — 0, in which case either two or all three are repeated, 

(ii) (no-Zeno region) A{rf) > 1 real root and 2 complex conjugate roots . 

The reason for the names for the two regions, Zeno and no-Zeno, will become clear when we 
examine the mixing behavior of random walks with dynamics prescribed by a dynamical matrix A 
that lies in one of the above regions. First, let us define the border between the Zeno and no-Zeno 
regions in terms of the parameters r x ,^f and A. This border is defined by the values where the 
discriminant equals zero: 



54A 2 



± 



54 7 A 2 



(48) 



When 7 < V12A the equation A(rf) = has no real solutions, and therefore we cannot have 
repeated roots in this parameter range. On the other hand, when 7 > \/l2A we have two real values 
for r\ that define the upper and lower boundary of the no- Zeno parameter region. Since r\ = 1 — r 2 
we get the following expressions: 



for V12A < 7 < 4A, repeated eigenvalues when r x = ± 



\ 



54A 2 



± 



12A 2 



54 7 A 2 



(49) 



for 7 > 4A, repeated eigenvalues when 



± 



\ 



54A 2 



( 7 2 - 12A 2 ) S 
54 7 A 2 



(50) 



Mathematically, the change in behavior at 4A exists because the minus branch of the expression 
under the square root in equation (49 1 becomes imaginary. The phase diagram in Fig. [I] shows 
the eigenvalue regimes as a function of the parameter ratio 7/A. Also, in this figure we explicitly 
show the regions in "r'-space corresponding to the two eigenvalue regimes for two specific values of 
7/A : y^IE and a/20- 

From this information about the eigenvalue phase diagram, we can make the following conclusions 
about the mixing behavior for various parameter regimes: 
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0.8- 
0.6- 



0.4- 




(a) Phase diagram for eigenvalues of the dynamical (b) r-space diagram of 

matrix A. The red (shaded) region indicates where eigenvalue types for matrix A 
eigenvalues are purely real, or what is referred to as when 7/ A = \/l5. 

the Zeno-region in the main text. 



(c) r-space diagram of 
eigenvalue types for matrix A 
when 7/ A = V20. 



FIG. 1: Eigenvalue regions of the matrix A in parameter space. 



(i) for 7 < Vl2A we have weak decoherence, so that for all decoherence axes (all r ) we have 
a pair of complex conjugate eigenvalues and therefore classical mixing and the possibility of 
instantaneous mixing. 



(ii) for 7 > V12A and the decoherence projection direction r x is such that \r x \ < 



V( 7 2-12A3) 3 



wc have two subcases: 



3 54A 2 ~ 54 7 A 2 

(a) if \/l2A < 7 < 4A and decoherence projection direction r x also such that \r x \ > 



I — 5452 — ^ 7 547A3 ~^~ tnen an eigenvalues are real, possibly repeated. There is no os- 
cillatory behavior in (a z )(t), and therefore no instantaneous mixing time exists. However, 
the walk of course has a classical mixing time as (er 2 ) t ~ , °° > 0. 

(b) if 7 > 4A all eigenvalues are real, possibly repeated. There is no oscillatory behavior in 
(a z )(t), and therefore no instantaneous mixing time exists. However, the walk of course 

has a classical mixing time as (a z ) t ^°°> 0. 

(Hi) for 7 > y/l2A and decoherence projection direction r x is such that it does not satisfy the 
conditions from (ii) and (ii-a) then we have a pair of complex conjugate eigenvalues and 
therefore classical mixing, and the possibility of instantaneous mixing. 

The fact that no finite instantaneous mixing time exists when all eigenvalues are real (for their 
restricted decoherence model) was interpreted by Alagic and Russell [2T] as an analogue of the 
Zeno effect where the quantum evolution of a system is hindered by its strong interaction with an 
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environment [3H]. Following this, we refer to the region where no instantaneous mixing time exists 
- i.e. regime (ii) above where all eigenvalues of A are real - as the Zeno-region, and the remainder 
of (7, A, r ) parameter space as the no-Zeno region. 

The presence of imaginary eigenvalues does not guarantee an instantaneous mixing time because 
it could still be possible that the timescale of the damping effect of the decoherence on (a z (t)) can 
be much faster the the characteristic period of the oscillation. We can formulate implicit conditions 



for the presence of instantaneous mixing by examining the expression from Eq. (38) for (a z (t)) when 
A has imaginary eigenvalues: 

(a z (t)) = A ie Xlt + A\e x '^ + A 3 e^ (51) 
where we have assumed Ai and A2 — A| arc the imaginary eigenvalues. Ai are the time independent 



coefficients of the exponentials in Eq. ( 38 ) . A necesarry and sufficient condition for instantaneous 
mixing is: 3 t s.t. (cr z (t)) < 0. Writing - ^! = | Ai le 1,6 * and Ai = \r.i + iXi A — ^ and using the fact 
that A3 = — Ai — \\ — 27, this condition can be written as: 

3t s.t. 2|Ai|e 3A " 1 *cos(6' + X IA t) < -A 3 , (52) 

which reduces down to: 

3t s.t. |^i|e 3AH lt cos(6»+ V3X R ,it) < -2(3A| (1 - Ad 2 + r 2 ) (53) 

where A 1 = |(4d 2 + 9X 2 R1 - r 2 z ) + i^(id 2 - 3X\ X - r 2 z (l + 4Xr,i)) and d = ^. For any value of the 
physical parameters, (7, A, r ), this condition can be checked for the existence of an instantaneous 
mixing time. 



C. Numerical simulations 



In this section we numerically evaluate (a z )(t) and calculate the classical mixing time for several 
parameter values and e = 0.001. These simulations are summarized in Fig. [2] which show mixing 
time for all values of ~r and three values of the physical parameter ratio 7/ A. We make several 
observations from these plots: 

• The mixing time can vary considerably as ~r is varied. Although the change in mixing time is 
generally smooth with changes in r , there are regions where the change is abrupt (e.g. around 
r x = ±1 when 7/ A = 1). This implies that the mixing time can potentially change drastically 
with the exact value of r . Hence it is very important to characterize the decoherence process 
accurately for determining mixing properties. A similar conclusion was arrived at by Strauch 
in Ref. [3] where he demonstrates that mixing behavior differs greatly depending on the 
decoherence model chosen. 

• The range of mixing times on the ~r*-sphere is smallest when 7/ A w 1. The range of mixing 
time diverges when this parameter ratio is very large or very small. This suggests an optimal 
7/ A parameter ratio where the interplay between Hamiltonian dynamics and decoherence is 
such that decoherence in any direction yields small mixing times. We will investigate this more 
thoroughly in the next section. 

• In the Hamiltonian dominated regime, where 7/ A < 1, we have fast mixing near the r x axis 
but long mixing times in the r y — r z plane. 

• In the decoherence dominated regime, where 7/ A > 1, we have fairly similar mixing times 
across nearly all values of ~r except for a small region around r = (±1, 0, 0). The size of this 
region shrinks as 7/ A increases, but the value of the mixing time in this region grows with 
the same parameter. However, note that when the decoherence is exactly along the r x axis we 
have very short mixing times as evident from the exact solution given by Eq. (refeq: :Xbit-dep) 
for this case. 
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(c) 7 /A = 5 

FIG. 2: Variation of single qubit mixing times with decoherence axis for several values of 7/A. 
The color at each point on the Bloch sphere corresponds to the mixing time for single qubit 
dynamics when the decoherence vector is {r x , r y ,r z ). The mixing times are in units of A and are 

calculated for e = 0.001. Note that (a z )(t) is invariant under negation of any coordinate of the 
decoherence axis (r x — > —r x ,r y — > —T y ,r z — > — r z ), and so the portions of the sphere that cannot 

be seen can be inferred. 



D. Optimal decoherence rate 



The numerical simulations presented in the last section suggested that the smallest mixing time is 
achieved for a non-zero value of 7/A. The simulations also suggested that the optimal 7/A is around 
the critical ratio 7/A = 1. These observations match with Kcndon and Trcgcnna's conclusions in 
Ref. |16j where they showed that some amount of decoherence can lead to faster mixing of quantum 
walks on a line and cycle, and faster hitting times of the quantum walk on a hypercube. Here we 
fully characterize the scaling of mixing time with decoherence rate (for the hypercube quantum walk) 
by numerically evaluating the mixing time for several fixed decoherence axes. Figure [3] shows how 
the mixing time varies with 7/A for several choices of ~r . In general, the curves are similar for any 
latitude in l^-space, that is, for any fixed 9. For a given 9 the mixing time versus 7/A curves (for 
various <j)) show maximum variation when 9 = n/2 (i.e. when ~r is in the x — y plane). Therefore 
we have shown these curves on a separate plot in Fig. [3] 

The notion of an optimal ratio 7/A is accurate for nearly all decoherence axes. And for nearly 
all decoherence axes, this optimal value is in the range 1 < 7/A < 5. However, when ~r is in the 
x — y plane, there is no finite optimal value for 7/A; the mixing time decreases continuously as 
7/A is increased for r in the x — y plane. We can gain intuition about this result by viewing the 
decoherence as a localizing phenomenon - it tends to localize the qubit state along the axis (on its 
Bloch sphere) defined by ~r vector. And the larger 7/A is, the faster this localization happens. 
For a qubit localized in the x — y plane, (er 2 ) = and hence a fast localization to this plane yields 
fast mixing. T he v ariation of the mixing time when ~r is within the x — y plane is an interesting 



feature of Fig. 3b As r approaches the x-axis (</? = 0) the mixing time evolution becomes closer 
and closer to exponential decay (with 7/A). However, away from the x-axis, the curves still show a 
local minimum around 7/A ss 1, but the global minima are still for 7/A — > 00. 
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— e=it/2, <j)=o 
e=n/2, 4>=jt/5 

--6=ji/2, <j)=3n/10 

- e=ji/2, (ji=jt/2 




10 20 30 40 50 

y/A 

(b) Characteristic behavior of mixing time versus 7/A for 
9 = tt/2. 

FIG. 3: Scaling of (log of) single qubit mixing time with the physical parameter ratio 7/A. The 
mixing time curves are shown for several values of 9 and ip, the two angular parameters of the 
decoherence axis. Figure (a) shows mixing time curves for < 9 < ir/2. We primarily only show 
curves for ip = in this parameter range for 9 because the behavior of mixing time versus 7/A for 
other values of ip is very similar: diverging mixing time as 7/A — > and 7/A — > 00. As 9 —* n/2 
these curves show more variation with ip, but they maintain this general shape. Figure (b) shows 
mixing time curves for 9 = tt/2 and several values of ip. 



VI. CONCLUSION 



We complete the picture that exists in the literature for quantum random walks on the hypercube 
with decoherence under the subspace projection decoherence model. This model of decoherence is 
the most physically realistic form when the quantum random walk is implemented using a register 
of qubits. 

The following are the important points of our work: 

• The instantaneous mixing time property of the quantum random walk does not necessarily 
disappear as the decoherence strength is increased. More precisely, Zeno dy nami cs prevails only 



in the precisely specified regions of the r , 7, A parameter space defined in VB Consequently. 



the hypercube quantum walk with decoherence does not always limit to its classical version as 
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the decoherence is increased. Depending on the direction of the decoherence vector it is possible 
to have oscillatory dynamics and instantaneous mixing persist as decoherence is increased. 

We showed numerically that, for almost all decoherence directions, a finite decoherence rate 
exists for optimally fast mixing. This optimal rate is approximately the same as the rate of 
Hamiltonian evolution. However, we also showed that for certain decoherence directions ( r in 
the x — y plane) no finite optimal decoherence rate exists and the mixing rate increases without 
bound as the decoherence rate is increased. This result is particularly relevant given recent 
results on quantum walk based modeling of excitation transport in biomoleculcs [39, 4TJ1I3T]. In 
these works it is argued that decoherence can lead to faster hitting times and walker diffusion, 
and this is explicitly confirmed for a simplified model of dephasing of the quantum walk. Our 
results suggest that the exact model of the decoherence matters greatly, and therefore, that an 
accurate model of the environmental interactions is essential to assess the merits or demerits 
of decoherence to excitation transport. 



Furthermore, we showed in section V A that randomizing the decoherence axis yields the 



depolarizing channel which exhibits both instantaneous and classical mixing regardless of the 
relative decoherence rate ^ . Yet another example of a quantum walk that does not decohere 
to a classical walk as the decoherence rate is increased. Therefore, introducing randomized 
decoherence may be an avenue for controlling the mixing behavior of the hypercube random 
walk. 



As detailed above, when implemented using a register of qubits, the hypercube random walk 
possesses an important symmetry that allows one to analyze its properties by simply considering the 
dynamics of a single qubit. Exploiting this property was crucial to deriving our results. Examining 
decoherence models on general quantum walks, especially ones without a lot of symmetry, is a 
difficult task. While we do not expect all our techniques to directly port to such situations of 
greater complexity, we do believe that the approach of converting the dynamics of a noisy quantum 
random walk into an equivalent continuous time dynamical map will be a fruitful one. 
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APPENDIX A: SINGLE QUBIT MASTER EQUATION 



The master equation, (16 1, is hard to solve in general but in our case the system Hamiltonian 
and decoherence operators are a sum of the tensor products that have a special structure. Each 
summand is a tensor product of elements only one of which is not the identity. We now show formally 
that this allows one to consider a combination of single qubit evolution equations. This calculation 
generalizes the calculation in |21j to deal with a decoherence axis that is along an arbitrary direction 
on the Bloch sphere. 

Vectorization is a technique that transforms any n x n matrix into a n 2 dimensional vector by 
stacking the transposed rows of the matrix on the top of each other. We will denote a vectorized 
matrix X as X v . A useful identity we will utilize involves the vectorization of the matrix product 
AXB: (AXB) V = (B T <g> A)X V . The action of unitary evolution on a density operator is UtpU}, 

and that is vectorized as (UtpU}) v = (Z7 t * ® Vt)p v = StP v , where St is the matrix form of the unitary 
evolution superoperator. Using this formalism it is straight forward to derive the vectorized picture 
of the master equation. Consider the discretized evolution given by Eq. ^ with r — * dt, and the 



qubit projection POVMs given by Eq. ( |14| 

Bt +dt = (1 - jdt)U dtQt Ul t + ydtJ2 E M k a [U dtQt ul t ]M k a (Al) 



fe a 



This can vectorized as: 



Ql +At = {l- 1 dt){U dtQt u\ t y + jdtV v (U dt gtUir (A2) 
= (1 - 7 *) (U* dt ® U dt )g v t + !dtV v (U* dt ® U dt )g v t (A3) 
= [(1 - 1 dt)+ 7 dtV v ][US t ®U dt ] Q v t (A4) 



where V v is the operator: 



l 



fe=l a=0 
n 

= J2 I ® 2(k ~ 1) ® pP5(^)®Po(^) + PJ(^)®iPi(^)] < 

k=l 

This transition equation for g defines the system dynamics at all times. Now, let g v t = S t g V Q where 



St is the propagator matrix for the dynamics. From Eq. (A4.I we know that g v t+dt = [(1 — jdt) + 
r ydt'P v ][U d ' t ® Udt]StQo, and since this is true for any initial state Qq, we get: 

St+dt = [(i-"fdt) + 1 dtr v ][u* dt ®u dt ]St 

= [(1 - jdt)l + "fdtT v ][(l+iHdt) <g> (l-iHdt)]S t (A5) 
= [l®l+idt(H®l-l®H)--fdtl®l+jdtV v ]S t , (A6) 

where we have expanded U d t — e~ lHdt to first order, and 1 = l® n . Taking into account that 
S t = ^ = lim di ^o St+d d t i T St we get the differential form: 

S t = [i[H ®l-l®H]-7l®l + IP ] S t (A7) 
Using H = A Yl=i l® {k - 1] ®a x ® we expand this as: 



iA[<7 x (g> I - I <g) <J X ] - 7I (g) 1+ \ ^ Tr «,2(n-l-fc) 

(r ) «) PoC? ) + PJ( "r ) ® Pi("r)] 



5 t (A8) 

AS t . (A9) 
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Since the initial condition is So = I® n ., the solution to this differential equation is St — e At : 



S t = exp 

n 

_ y^®2(fc-i) 



Emg>2(fe-1) a f / iAka, ® I - I ® CrJ - 7I ® 1+ 
l 7Po( r ) ® PoC? ) + ("r ) <8 Pifr )] 

_k— 1 v 



1 expt 



fc=i 



iA]^ g) I - I ® crj - 7I <g> 1+ 
7pP5( r ) <g> P„("?) + PJC?) <g> Pi("r )] 



i-i-fc) 



ri-l-fc) 



(A10) 



= [expt {iA^ ® I - 1 <g> a x ] - 71 1 + 7pP?("r ) <8> P ("r ) + ® Pi( "r )] }]® n 

= [S t ] 8n 

Therefore dynamics of the system is tensor product of individual qubit dynamics St — e -4 * . The 
single qubit generator can be simplified as follows: 



= iA[a x ® I - I ® a x ] - 7I ® I 



+ I[(l + "r ."ct)* ® (1+ ■ ~a) + (l - r* • ~a)* <g> (l 



= iA[<7;c 



7, 



7- 



r • er (g) r • a . 



(All) 

(A12) 
(A13) 



It can be easily confirmed that this is the generator for the single qubit dynamics described by Eq. 
(17 1 once it has been vectorized. 



APPENDIX B: DEPOLARIZING CHANNEL THROUGH THE RANDOMIZED 

DECOHERENCE AXIS 

The generator for the single qubit dinamics when the decoherence axis is randomized is: 



A d = i&{a x ®l-l®a x )~ 1 l®l+ & [P5(r)®P (r ) + P*(r )(g)Pi(r)]ds . (Bl) 

47T J S 2 

We can carry out this integral to get the following: 



A d = 



iA(a x <g) I - I ® cr x ) - 7I (g> I 



dr 



iA(a x ® I — I ® a x ) — 7I ® I 

+ 8? lej 1 ® 1 + ^ ■ ^* ® • -?)tf? 

= iA^ (8 I - I ® 0- x ) - ^1 <g> I + ^((7 K <8S CT^ + <7* «) <T y + <J Z ® <T Z ). 

where we have used r x — sin 9 cos 95, r. y = sin#sin<ys, r z = cosO. This leads to the solution for the 
single qubit dynamics: 



P v (*) = e**Y (0) 



(B2) 



Changing the basis for Ad to the eigenbasis basis, exponentiating, and returning to the original 
basis we get: 



Pit) 



if l + e^cos(2Ai) -ie^^r sin(2Ai) \ 
2l ie"^ sin(2Ai) 1 - cos(2Ai) J 



(B3) 



with ,o(0) = |0)(0|. The value of p(t) QQ and /9(t) n determines the probability of measurement in 
basis {|0), |1)}. The eigenvalues of our operator a(A d ) = {0,-^,-^(7 - 3i), ^(7 + 3i)} 
determine the probability distribution. The expressions for p(t) Q0 and p{t) xl show that regardless 
of the rate of decoherence instantaneous mixing exists. 



